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Abstract 

In this present work, we try to build up a cosmological model using a non-canonical 
scalar field within the framework of a spatially flat FRW space-time. In this context, we 
have considered four different parametrizations of the equation of state parameter of the 
non-canonical scalar field. Under this scenario, analytical solutions for various cosmological 
parameters have been found out. It has been found that the deceleration parameter shows 
a smooth transition from a positive value to some negative value which indicates that the 
universe was undergoing an early deceleration followed by late time acceleration which is 
essential for the structure formation of the universe. With these four parametrizations, the 
future evolution of the models are also discussed. It has been found that one of the models 
(Generalized Chaplygin gas model, GCG) mimics the concordance ACDM in the near future, 
whereas two other models (GPL and JBP) diverge due to future singularity. Finally, we have 
studied these theoretical models with the latest datasets from SN la + H{z) + BAO/CMB. 

keywords: Dark energy, Non-canonical scalar, Parametrization, Data analysis 


1 Introduction 

Recent cosmological observations [1, 2, 3, 4, 5, 6] strongly suggest that our universe is presently 
accelerating. In literature, there has been a number of theoretical models to explain the ori¬ 
gin of this acceleration mechanism. In this context, the most accepted idea is that an exotic 
component of the matter sector with large negative pressure, dubbed as “dark energy” (DE), is 
responsible for this accelerated expansion of the universe. DE also makes up about 73% of the 
total energy budget of the universe at present epoch. However, understanding the origin and 
nature of DE is still a challenging problem in modern cosmology. A number of models have been 
proposed phenomenologically as DE models, such as quintessence (canonical scalar held) [7, 8, 9], 
phantom [10, 11, 12], k-essence [13, 14, 15], Chaplygin gas [16, 17, 18], /(R)-gravity models 
[19, 20, 21, 22, 23] and so on. The simplest theoretical candidate of DE is the vacuum energy 
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with a constant equation of state (EoS) parameter uj = —1, but it suffers from cosmological con¬ 
stant problem [24, 25]. The dynamical nature of dark energy also introduces a new cosmological 
problem, namely, “coincidence” problem [26]. One alternative to the coincidence problem are 
coupled dark energy models where DE interchanges energy with the dark matter (DM) by means 
of a coupling term [27, 28, 29, 30, 31, 32]. Though a number of theoretical models have been 
proposed, none of them provides a satisfactory solution to all the problems. Hence, there is still 
a need of an appropriate model to explain current observations. 

In the proposed models of dynamical dark energy, the EoS parameter is usually considered 
to be evolving with time. In this context, a large number of parametrizations of DE equation of 
state have been proposed [33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44], which could provide solu¬ 
tions to a number of cosmological problems. However, most of these analysis have been carried 
out for canonical scalar held models of DE. Recently, non-canonical scalar held models are also 
being studied as a candidate for DE. In the non-canonical scalar held models, the kinetic part of 
the scalar held is modihed and it has been found that an accelerated expansion can be achieved 
through these modihcations. Originally, Armendariz-Picon et al. [45, 46] proposed this scenario 
to explain inhation at high energies where the non-canonical scalar held efficiently plays the role 
of inflaton (for review on this topic, see [45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59]). 
Later, Chiba et al. [13] introduced this scenario for dark energy models. In our recent work [60], 
we have studied an interacting non-canonical scalar held model with a constant EoS parameter 
for the scalar held. The model was however restricted in the sense that a constant EoS parameter 
does not provide a general framework. Motivated by the above facts, in this present work, we wish 
to consider a varying EoS parameter for the non-canonical held. We consider four popular DE 
parametrizations for a non-canonical scalar held model in order to explain the late-time scenario 
of the universe; the parametrizations considered are Chevallier-Polarski-Linder parametrization 
[42, 43], Jassal-Bagla-Padmanabhan parametrization [44], Barboza-Alcaniz parametrization [38] 
and Generalized Chaplygin Gas parametrization [16, 17, 18]. 

The features of these various parametrizations have been discussed in details in the next sec¬ 
tion. We have obtained the expressions for diherent relevant cosmological parameters, such as the 
deceleration parameter, density parameters of the scalar held and matter held for each model, 
and have shown that it is possible to have late time accelerated expansion of the universe for 
each of these choices. We have also compared the results with standard canonical scalar held 
models considering these parametrizations. Furthermore, in this paper, we have discussed about 
the future evolution of the universe for these EoS parametrizations and we have found that one 
of these parametrizations (GCG) behaves as standard ACDM in far future where as GPL and 
JBP models fail to provide information about the far future and are valid till z > —1. BA model 
however can provide information regarding the entire evolution history of the universe till z = —1. 
Finally, we have studied the constraint on the EoS parameter using the combination of SN la -|- 
H{z) + BAO/CMB dataset. 

The paper is organized as follows. In section 2, we have described the basic theoretical frame¬ 
work for the non-canonical scalar held model of a hat FRW universe. We have then solved the 
governing dynamical equations for this toy model using four diherent types of DE parametriza¬ 
tions of the EoS parameter. It has been found that the resulting cosmological scenarios are in 
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good agreement with the current observations in each case. In section 3, we have obtained the 
observational constraints on this model parameters using SN la + H{z) + BAO/CMB dataset. 
Finally, some conclusions are presented in the last section. 


2 Field equations and their solutions 


The general action for a scalar held model (with SttG = c = 1) is given by 


A 


yj—gd'^x 


— + £(0, A) 


+ Sr, 


( 1 ) 


where R is the Ricci scalar curvature, C{(j)^X) is the Lagrangian density which is an arbitrary 
function of the scalar held 0 and its kinetic term X. The kinetic term X is dehned as A = 
\d^(j)d^(l) = 10^ for a spatially homogeneous scalar held and Sm represents the action of the 
background matter held. 

Varying this action with respect to the metric gives the Einstein held equations as 

1 dC 

R,u - -g,uR = -^d,(t>dA - + T- (2) 

where represents the energy-momentum tensor of the matter held which is modeled in the 
form of an ideal perfect huid and is dehned as 


Tny (Pm T Pmg^p 


( 3 ) 


where pm is the energy density, pm is pressure of the matter held respectively and is the four- 
velocity of the huid. Secondly, variation of the action with respect to the scalar held (f) gives the 
equation of motion for (p as 


dx dx^ 


+ 


rrdC 

-h 

, dX 


d^C 


dC 


dXd(^ ^ d(j) ^ 


( 4 ) 


The energy density (p^) and the pressure (p^) of such a held is given by 


Pep 



2X-C, 


P<p = ^ 


( 5 ) 


In general, the Lagrangian density for a scalar held can be written as [61] 


C{<P,X) = f{<P)F{X)-V{<P) 


( 6 ) 


where V{(p) is a self-interacting potential for the scalar held (p, F{X) is an arbitrary function of 
X. When f{(p) = 1 and F{X) = X the Lagrangian (6) reduces to the quintessence Lagrangian. 
It describes k-essence when V{(p) = 0 and phantom scalar held when f{(p) = 1 and F{X) = —X. 
In case of the phantom held, the sign of its kinetic term “V” is opposite compared to the action 
for a canonical scalar held. 
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The Lagrangian density for a general non-canonical scalar field is given by [59] 


C{(t>,X) = F{X)-V{(t>) 


( 7 ) 


These type of scalar held models with non-canonical kinetic term have received huge attention 
recently. Unnikrishnan et ah [57] have showed that for such models, the slow-roll conditions can 
be more easily satished compared to the canonical case. The non-canonical scalar held models 
are also found to be able to generate inhation in the early epoch [45, 46, 47, 48, 49, 50, 51, 52, 
53, 54, 55, 56, 57, 58, 59]. These attractive features of a noncanonical scalar held motivated us 
to study the features of this model in the context of dark energy. In this present paper, we have 
considered a Lagrangian density of the following form 

C(4>,X) = X'-‘-V(4,), x = ly (8) 

which can be derived from the general form of Lagrangian density [57, 62, 63] 


C{4>,X) = X 



V(4>) 


( 9 ) 


for a' = 2 and Mpi = -?= = 1. 

V ottO 

It must be noted that the above equation reduces to the well known Lagrangian density for a 
canonical scalar held model when a' = 1. It is also worth mentioning that the Lagrangian (8) 
dihers from the Lagrangian of k-essence and phantom models in the sense that for phantom held 
kinetic energy term itself is negative and here the potential term is non-zero. This type of La¬ 
grangian has also been considered in our earlier work [60]. In the subsequent sections, we try to 
build up an accelerating model for the universe in which the non-canonical scalar held will play 
the role of dynamical dark energy. We are basically interested to study how the dynamics of the 
non-canonical dark energy model gets ahected by the various parametrizations of DE equation of 
state parameter. The behaviour of canonical scalar held models are very well studied for diherent 
parametrizations of DE EoS. However, the ehect of these parametrizations have not been studied 
for a non-canonical scalar held sector. 


The energy density and pressure associated with this Lagrangian density can be obtained from 
equations (5) and (8) as 

+ V {(p) (10) 

= -V{<P) ( 11 ) 

The metric for a homogeneous, isotropic and spatially hat ERW model of the universe is charac¬ 
terized by the following line element 

ds"^ = dt^ — a^{t)[dr‘^ + r'^dO'^ + r^sin^Odcj?] ( 12 ) 

where a(t) is the scale factor, normalized so that at present a(t)|t=to = 1 and t is the cosmic time. 
The Einstein held equations for the space-time given by equation (12) with matter in the form of 
pressureless perfect huid takes the form, 

= p,„ + + V(4,) 


(13) 
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2H + ?>H^ = --^^ + V{<t>) (14) 

/9^ + + p<^) = 0 (15) 

pm + ^Hpm = 0 (16) 

Here an overdot indicates differentiation with respect to the cosmic time t. Among the above 
four equations (equations (13)-(16)), only three are independent equations with four unknown 
parameters H, pm, 0 and V{(j)). So we still have freedom to choose one parameter to close the 
above system of equations. For the present work, we consider various functional forms of the EoS 
parameter for the scalar held. 

It is well known that the parametrization of DE equation of state plays an crucial role in 
understanding the nature of DE component. In general, the EoS parameter can be parametrized 
as, 

u^{z) = — = uo +uif{z) (17) 

p<f) 

where uo, ui are real numbers and f{z) is a function of redshift z. It may be noted that the 
standard hat ACDM model is represented by this parametrization with the choice of Wq = —1 
and ooi = 0. In fact, many functional forms of f{z) have been considered in literature [33, 34, 35, 
36, 37, 38, 39, 40, 41, 42, 43]. In this present work, we shall use four popular parametrizations of 
u^{z) to study the behavior of the deceleration parameter q{z) of this non-canonical scalar held 
model. 

2.1 Chevallier-Polarski-Linder (CPL) parametrization 

Among various parametrizations, the CPL parametrization (for details, see Refs. [42, 43]) is one 
of the most popular ones and is given by 


j^{z) = Wo + wi(l - a) = Wo + wi 


where 2 ; = - — 1 is the redshift, Wo represents the current value of u^pi^z) and the second term 
accounts for the variation of the EoS parameter with respect to redshift. In this present model, 
we have considered CPL parametrization of EoS parameter because this parametrization has the 
advantage of giving hnite in the entire range, 0 < 2 < cxd. 

The solution for p^ from equation (15) is obtained as 

P<t>{z) = P<t>oi^ + (19) 

where, cti = (1 + wq + wi) and P(^o is an integrating constant. From equation (16), we have the 
expression for energy density of matter as 


pm{z) = Pmo(l + zf (20) 

where, pmo is an integrating constant. From equation (13), the Hubble expansion rate can also 
be written as 

H\z) = Hi b™o(l + zf + D<^o(l + (21) 
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Here, Hq is the Hubble parameter at the present epoch, and are the density- 

parameters at the present epoch of the matter and scalar held respectively. 


The deceleration parameter is dehned as, q = — = — (1 + ■^)- 


following form 


,.13 
qiz) = - + - 


^0 + (i^J 


1 + k( 1 + z)(3-3«i)e( i+z ) 


For this model, q takes the 


( 22 ) 


where, k = = 7 ^. o < 0 indicates accelerated phase of the universe while o' > 0 indicates 

a decelerated phase of expansion. From hgure (la), we see that q decreases from positive to 
negative value for suitable choices of model parameters. 



Figure 1: a) Plot of q as a function of z (upper panel) and b) Plot of VL^ (dashed curve) and 
(solid curve) as a function of z (lower panel). This is for k = Wq = — 1 and Ui = 0.1. 

Here, ai = (1 + coq + oji). 


For this model, the evolution of the density parameters of the matter and scalar held are 
obtained respectively as. 


^m{,z) 


1 

1 + i(l + 2;)3"i“3g(- i+z ) 


(23) 


( 3ui,z \ (24) 

1 + k( 1 + 2;)3-3oigt i+z ) 

which further yields, Qm(z) + f2</,(z) = 1. Figure (lb) shows the plot of density parameters 
for the scalar and the matter held as a function of ;s. This graph shows that starts domi¬ 
nating over flm at around z ~ 0.4. This result is compatible with the observational results [64, 65]. 



However, the model presented here is restricted because u,p(z) diverges when z —>■ —1 i.e., 
this model cannot predict about the future evolution. So, this particular toy model is capable of 
describing the evolution history of the universe from the past to the near future upto z > —1 but 
can not predict about the evolution beyond that limit. 
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2.2 Jassal-Bagla-Padmanabhan (JBP) parametrization 

Recently, Jassal et al. [44] extended the above parametrization to a more general case: 

Uj^{z) = ujq + (25) 

For the present model, we choose p = 2. It must be noted that the EoS parameter o;^ ~ ojq at 
both high and low redshifts for p = 2. Also, one can obtain the widely used CPL parametrization 
of EoS from equation (25) for p = 1. For the JBP parametrization, using equation (15), the 
expression for can be obtained as 

/ 3ujiz^ \ 

p4z) = p^o(l + z)3d+-o)elFTTiFj (26) 


where p^o is an integrating constant and represents the present value of the scalar held density. 
The Hubble parameter for this model takes the following form 


= Hi 




f2mo(l + + f2</)o(l + 


(27) 


In this model, we express deceleration parameter q as 


q{z) 


1 

2 


3 

2 


1 + «:(! + ^)-3^oe"V2(i+-)V 


The corresponding density parameters are now given by 


( 28 ) 




Figure 2: a) Plot of q as a function of z (upper panel) and b) Plot of flm (dashed curve) and 
(solid curve) as a function of z (lower panel). Both the plots are for n = wq = — 1 

and LJi = 0.1. 


1 


/ 3ujiz^ \ 

1 + 1(1 + 




(29) 
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(30) 
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/ 3ujiz^ 

1 + «(! + ^)-3a;oe"V5aT^ 

Figure (2a) shows the plot of q(z) as a function of This plot clearly shows the transition of g 
from the decelerating to the accelerating regime at z ~ 0.8. The evolntions of Qm and against 
z are shown in fignre (2b). The plots are for wq = and wi = 0.1. In both the graphs, the 
resulting cosmological scenarios are in good agreement with observations. For the JBP model 
also, uj^(z) diverges as z —1 and thus future evolution can not be predicted. 

2.3 Barboza-Alcaniz (BA) parametrization 

The next parametrization considered in this paper was proposed by Barboza et al. [38], which 
has the following fnnctional form 



u;^(z) = ujo + czi 


z{l + z) 
1 + z'^ 


(31) 


where Ufi,{z = 0) = wq (the present valne of the EoS parameter), oji = ^1^=0 (which measnres 
the variation of the EoS parameter with z), u^{z = oo) = uj^ + ui and the EoS parameter rednces 
to u^{z) = Wo + at the low redshift ( 2 : << 1). It is remarkable that the BA parametrization 
does not diverge like CPL model when z ^ —1. 


In this model, Pij){z) becomes 

p*(^)=p«(i+.y“+"“'(i+zt^ 


Now the eqnation (13) can be written as 


H^z) 


Hi 


f ^ mo(l + z )^ + fl 0 o(l + ^)^*'^''~‘^°^(1 + z '^) 2 


(32) 


(33) 


In this case, the deceleration parameter q{z) can be expressed as 


q{z) 


1 

2 


3 r + 

2 1 + k{ 1 + z)~^‘^^ {1 + z‘^)~ ^ 


(34) 


Furthermore, one can express the density parameters of the matter and scalar held respectively 
as ^ 

f^m(2^) ) 3iirr 

1 + i(l + 2;)3 ‘^o(1 + 2:2)~ 

1 

3u;i 

1 + k {1 + -f- 2 ( 2 ) — 

Figure (3a) shows the evolution of the deceleration parameter with redshift z. It is evident from 
hgnre (3a) that the universe is presently nndergoing an accelerating phase of expansion (g < 0). 
Fignre (3b) shows that the density parameter VLrn increases with 2 , whereas decreases with 2 :. 
This featnres of q{z), and are consistent with the present day observations. 
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Figure 3: a) Plot of q against z (upper panel) and b) Plot of (dashed curve) and (solid 
curve) against z (lower panel). Both the plots are for k = cvq = —1 and Ui = 0.1. 


2.4 Generalized Chaplygin Gas (GCG) parametrization 

It is well known that the generalized chaplygin gas (see Refs. [16, 17, 18]) behaves like dark 
matter in the past and it behaves like cosmological constant at present. Motivated by this idea, 
in this paper, we are interested to describe the late-time dynamics of the universe produced by 
the GCG. For this purpose, we have assumed that the universe contains both the dark matter 
and the GCG. Additionally, we have also considered another interesting possibility where the 
non-canonical scalar held 0 plays the role of GCG to explore the late time cosmic scenarios. The 
GCG equation of state is described by [16, 17, 18] 

Pi. = -4 ( 37 ) 

P% 

where A is a positive constant and a is another constant in the range 0 < a < 1. The original 
chaplygin gas corresponds to the case a = 1 [16]. By inserting equation (37) into the energy 
conservation equation (15), one hnds that the density of the scalar held cj) evolves as 




A + R(1 + z)3(^+“) 


1 

ITT^ 


(38) 


where, B is an integration constant. Equation (38) can be re-written in the following form 


— P4>o As -|- (1 — As)(l -|- 


a) 


(l + oc) 


(39) 


where, for simplicity, we have dehned Ag = and = {A + B) (i+“) is the present value of 
the energy density of the GCG. To ensure the hnite and positive value of we need —1 < a < 1 
and 0 < As < 1. 


In this case, the Hubble parameter is given by 


= Hi 


Gmo(l + ^)^ + G(^o(As -1- (1 — As)(l -|- 2 ;)^d+“)^ (l+a 


1 H 

7 


( 40 ) 
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The corresponding expression for the EoS parameter is given by 

^ “^ + (l-^)(l + z)3(l+“) 

Like earlier mentioned three models, the EoS parameter of the GCG also depends on two inde¬ 
pendent model parameters and a) along with redshift 2 ;. At present epoch, the above EoS 
parameter becomes, Uip{z = 0) = —Ag. It is interesting to note that the GGG will behave like 
pure cosmological constant when we put A,, = 1. 


The deceleration parameter q can be written as 

_ A, 

A^ + (l-Aa)(l+^)3(i+^) 
k(1+2)3 

- 1 

[>1s-I-(1-As)(1+2)3(i+“)] (i+“) 

where k = ^ 

P<^0 s 

In this case the density parameters have the following form 





z z 


( 42 ) 


Figure 4: a) Plot of q as a function of z (upper panel) and b) Plot of flm (dashed curve) and 
(solid curve) as a function of z (lower panel). Both the plots are for Ag = 0.9, a = —0.5 and 

^ 0.27 
^ 0.73 ■ 






1 


1 + 

[A„+(1-A„)(1+2)3(i+“) 

1 

li+SJ 

k(1+2:)3 

1 


1 + 


-1- 


[A„ + (1-71„)(1+2)3(1+“)] (1+“) 


(43) 

(44) 


Figure (4a) shows the evolution of q{z) vs. redshift z for Ag = 0.9 and a = —0.5. In fact, at low 
redshift, the transition of q{z) from decelerating to accelerating regime depends upon the choice 
of the parameters Ag and a. Also, the evolutions of flm and fl^f) against 2 : are shown in figure (4b) 
for the earlier mentioned same chosen values of Ag and a. 
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2.5 Comparison between canonical and non-canonical scalar field mod¬ 
els for the above parametrizations: 

For all these models, the relevant potential for the scalar held (j) in terms of redshift z can be 
written as (from eqnations (10) and (11)) 


V{z) = -{I - 3u^{z))p^{z) 


which immediately gives 


Vcpl{z) — Vq (1 — 


3uJiZ 

1 + z 


o 3a;i z 


(45) 


(46) 


VjBp{z) 

Vba{z) = Vo 


: H (^1 - Swo - 

l-3uo- (1 + (l + z^) ^ 


(47) 

(48) 


Vgcg{z) = Vo{4As + (1 - d,)(l + X [a, + (1 - 7l,)(l + z)3d+' 


ha) 


(l + c) 


( 49 ) 


where Vq = _ Vcpl-, Vjbp, Vba and Vgcg are the potential for the CPL, JBP, BA and GCG 

models respectively. Here, Hq and f2(^o represent the present day valnes for the Hnbble parameter 
and the dark energy density parameter respectively. 

Adding eqnations (10), (11) and replacing 0 = aH^, one can obtain the general expression for 



<l> 


Figure 5: This figure shows the variation of the potential V with for non-canonical scalar field 
model by assuming wq = — 1, (Ui = 0.1 for CPL (thick curve), JBP (dashed curve) and BA (dotted 
curve) parametrizations and Ag = 0.9, a = —0.5 for GCG (thin curve) parametrization. All the 
plots are for the parameter choices Vt^o = 0.73, 0o = 0.1 and Ho = 72 kms~^Mpc~^ 

the non-canonical scalar held as, 

^ [{l + U^{z'))p^{z')]~^ / 

{l + z')H{z') 


(f){z) = (j)o + 


0 


( 50 ) 
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where 00 is an arbitrary integration constant. 


Now, we will focus on the extensively studied canonical scalar held case, in which the La- 
grangian density is obtained from equation (7) as 

C{<P,ano,X)=X-V{(P) (51) 

The energy density and pressure for the scalar held {(jicano) are given by 


Pcano Vcano{(pcano) i Pcano ^^cano Vcano{(pcano) 

where, Vcano is the potential of the canonical scalar held. The EoS parameter Ucm 
dynamical variable which gives the continuity equation (15) in an integrated form 


(52) 


_ Pcano 

pcano 


IS a 


pcano — Po^^P 


1 + ‘^'canoi^') 

1 + z' 


dz' 


where po is an integration constant. The Friedmann equation then becomes 

n 3 , /I 0 5 1 +'^cano(^0 


HLnoiz) = + ^)^ + (1 


In this case, the expressions for Vcano and (f)cano can be written as 



1 + ^' 


dz' 


hcano('2') 2^^ ^cano(^^')') Pcano(^^') 


and 


(z) — 00 + 


[(1 + ^cano {z'))p canoi,^ )] ■ 


■dz' 


(53) 


(54) 


(55) 


(56) 


(1 + z')H^no[Z) 

In the previous subsection, we have discussed various parametrizations of the EoS parameter for 
the non-canonical scalar held. Now, we will consider these parametrizations for canonical scalar 
held models to compare their behavior with the non-canonical scalar held models. It deserves 
mention that one obtains the same expressions of Pcano{z) and Hcano{z) for both canonical and 
non-canonical scalar held models. But, the expressions for potential associated with the scalar 
held 0 will be diherent for canonical and non-canonical scalar held models (see equations (45), 
(50), (55) and (56)). 

The expressions for the potential V{z) and (j){z) are very complicated and it is very difficult to 
express V in terms of 0. So, we have solved equations (50) and (56) numerically and have plotted 
E as a function of 0 (see Figures 5 and 6). In hgure 6, for each panel, the dashed curve shows the 
evolution of the potential Vcano{4>cano) for each parametrization whereas the solid lines represents 
the evolution for corresponding non-canonical case. It has been found that for canonical case, the 
slope of the potential is quite hat and (pcano is almost constant {(plano << Kano) throughout the 
evolution, which yields pcano ~ Vcano = constant, however, the case is diherent for a non-canonical 
scalar held model. Figure 6 also shows for each panel, that for non-canonical case, the trajectory 
of the potential V (0) changes very slowly at early epoch, but it starts increasing with 0 in such 
way that the potential term dominates over the non-canonical kinetic term (0"^) at late times 
independent of initial conditions and thus provides acceleration. So, the non-canonical scalar 
held exhibits an interesting property of the potential, which can provide a possible solution to 
the coincidence problem. This new feature occurs due to the non-canonical kinetic term present 
in the Lagrangian (8). 
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(f) (f) 




<l> 


0 


Figure 6: This figure shows the variation of the potential V with by assuming uq = —1, oji = 0.1 
for CPL, JBP and BA parametrizations and Ag = 0.9, a = —0.5 for GCG parametrization. The 
dashed curve represents the trajectory of the potentials for the canonical scalar field, as shown in 
each panel. All the plots are for the parameter choices = 0.73, Vtmo = 0.27, fio = 0.1 and 
Ho = 72 kms~^Mpc~^. 


3 Observational Constraints on 

In this section, we shall fit these four parametrized theoretical models with recent observational 
datasets, namely Type la Supernovae (SN la), measurements of Hubble parameter (H(z)), bary- 
onic acoustic oscillations (BAO) and cosmic microwave background (CMB) observations. Al¬ 
though the procedure for calculation of individual function is quite well-known, we briefly 
mention the same for completeness. 


The total joint data analysis is dehned as 


2 _ 2 I 2 I 2 

Xtotal — XSN ' Xh ' XbAO/CMB 


(57) 


where the individual yfi for each dataset is evaluated as follows. 

First, we have used the latest observational dataset of SN la [66] of 580 data points. In order 
to calculate xfi^ for SN la data we follow the procedure described in Refs. [67, 68, 69, 70]. For 
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SN la dataset, the function is constructed as 


2 . 

XSN — ^ ~Q 


(58) 


where A, B and C are dehned as follows 


i=l 

® = -3- 


(59) 

(60) 


and 


580 1 

i=l 


(61) 


where represent the observed and theoretical distance modulus respectively and Uj rep¬ 

resent the uncertainty in the distance modulus. 


Next, we have used the 29 data points of H{z) [71, 72, 73, 74, 75, 76, 77, 78, 79] in the redshift 
range 0.07 < 2 ; < 2.34. These values are presented in table (1). To complete the dataset, we fixed 
the valne of Hq from [80]. For this dataset, the function is dehned as 


2 ^[h°^^{Zi)-h^\zi)f 

^ 0-2 (-y\ 


(62) 


where h = is the normalized Hnbble parameter and an is the error associated with each data 
point. In above eqnation snbscript “o5s” refers to observational quantities and subscript “th” 
refers to the corresponding theoretical ones. 


We have also used BAO [81, 82, 83] and CMB [84] measnrements dataset to obtain the 
BAO/CMB constraints on the model parameters. In order to calculate Xbao/cmb function, we 
follow the procednre described in Ref. [85]. We use the measurement of “acoustic scale” I a 
provided by CMB, which is dehned as 


Ia = tt 


dAjz^) 

rs{z^) 


(63) 


where dA^zA) = JA is the comoving angular diameter distance, rs{zA is the comoving sound 
horizon at the photon-decoupling epoch and ~ 1091 is the decoupling time. In this work, we 
have used I a = 302.44 ± 0.80 [84]. For this analysis, we have used six data points from 6dFGS 
[81], SDSS LRG [82] and WiggelZ [83] surveys. Following Ref. [85] and combining these results 
with [84], we obtained the BAO/CMB constraints on the model parameters (see table 2). 

The for the BAO/CMB data is given by 

Xbao/cmb = X C X (64) 
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where, X (transformation matrix) and C ^ (inverse covariance matrix) which are formed using 
different functional form of For details refer to Goistri et ah [85]. 

The GPL, JBP and BA model have three free parameters, namely, Gmo (or = 1 — G^o), <^o 
and Ui- In this case, the conhdence region ellipses in the uq — ooi parameter space can be drawn by 
hxing VLrno to some constant value. So, we have done analysis by hxing Gmo (the present value 
of the density parameter of the matter held) to 0.26, 0.27 and 0.28 for those dataset. Hence, we 
can now deal with only two free parameters (wq, i^i) and will obtain the observational bounds on 
this parameter from analysis of the combined dataset (SNIa + H(z) + BAO/GMB). With this 
we have plotted la (68.3%) and 2a- (95.4%) conhdence contours on uq — Ui parameter space for 
various DE parametrizations. Similarly, the GGG model has three free parameters, namely, GmO) 
As and a. In this case, we have plotted la and 2a conhdence contours on Ag — Gmo parameter 
space by considering a = —0.5. The best-ht values of the parameters are obtained by minimizing 
y^. In hgure (7) and (8), the large dots represent the best ht values of the model parameters 
and the small dots represent the chosen values of these parameters in our analytical models (as 
mentioned in previous section). In our analytical models, we have chosen wq = = 0.1 (for 

GPL, JBP and BA model) and Ag = 0.9, Gmo = 0.27 (for GGG model) as we are interested to 
study and understand the ehect of individual DE parametrizations for some hxed values of model 
parameters. For this combined datasets, the chosen values of Uq and Ui (A^ and a) are found to 
be well within the la conhdence contour. The observational bound on these model parameters 
as well as their best ht values are presented in table 3 & 4. For each model, we notice from table 
3 & 4 that the best ht value of the present EoS parameter is very close to —1 which is consistent 
with the recent observations [87, 88]. The standard hat AGDM model (wo = —1 and oji = 0) 
corresponds to the intersection point of the dashed lines as plotted in hgure (7), and it is evident 
from the hgures that the AGDM model is always inside the la conhdence contour for all these 
parameterizations. 

In hgure (9), we have shown the behavior of the EoS parameter u!^(z) = , for four 

diherent parametrizations using the best ht values of the model parameters (as listed in Table 3 
& 4) for this dataset. It may be noted that the forms of EoS parameter (or potential) are diherent 
for diherent models. But, uj^(z) ~ —1 at present epoch (z = 0 or a = 1) for GPL, JBP & BA 
parametrization model, whereas for GGG model, uj^(z) is around —0.9 at present but approach 
— 1 value as evident from hgure (9). Gonsequently, this feature suggests that (fA « V{(f>) i.e., 
the potential term dominates over the non-canonical kinetic term to accelerate the cosmic expan¬ 
sion at present. It is also evident from hgure (9) that the non-canonical scalar held behaves like 
phantom dark energy < —1, [10, 11, 12]) at present epoch for GPL, JBP and BA models. 
However, it deserves mention again that these plots of u^{z) are for the best-ht values of the 
model parameters and the results may change for the combination of other datasets. In this 
context, we would like to mention that the hypothetical phantom models meet several difficulties 
such as classical and quantum instabilities due to its negative kinetic energy and momentum as 
studied in refs. [89, 90, 91, 92]. In spite of that the recent observational data allow the possibility 
of the phantom EoS in the near past or in the near future [44, 93, 94, 95, 96, 97]. So, from 
observational viewpoint, the phantom held one of the possible candidates for dark energy and 
can not be ruled out completely. Interestingly, it has been found that for the BA model, the EoS 
parameter does not deviate much from —1 in future (upto 2 : = —1) and provides information 
regarding the complete evolution history of the universe but for GPL & JBP models however, the 
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Table 1: H{z) measurements (in unit [km s ^ Mpc ^]) in the redshift range 0.07 < 2 : < 2.34 and 
their errors._ 


z 

H{z) 

CTh 

Reference 


69 

19.6 

Zhang et ah [71] 

0.1 

69 

12 

Simon et ah [72] 

0.12 

68.6 

26.2 

Zhang et ah [71] 

0.17 

83 

8 

Simon et ah [72] 

0.179 

75 

4 

Moresco et ah [73] 

0.199 

75 

5 

Moresco et ah [73] 

0.2 

72.9 

29.6 

Zhang et ah [71] 

0.27 

77 

14 

Simon et ah [72] 

0.28 

88.8 

36.6 

Zhang et ah [71] 

0.35 

82.7 

8.4 

Chuang et ah [74] 

0.352 

83 

14 

Moresco et ah [73] 

0.4 

95 

17 

Simon et ah [72] 

0.44 

82.6 

7.8 

Blake et al. [75] 

0.48 

97 

62 

Stern et al. [76] 

0.57 

92.9 

7.8 

Samushia et al. [77] 

0.593 

104 

13 

Moresco et al. [73] 


87.9 

6.1 

Blake et al. [75] 

0.68 

92 

8 

Moresco et al. [73] 

0.73 

97.3 

7 

Blake et al. [75] 

0.781 

105 

12 

Moresco et al. [73] 

0.875 

125 

17 

Moresco et al. [73] 

0.88 

90 

40 

Stern et al. [76] 

0.9 

117 

23 

Simon et al. [72] 


154 

20 

Moresco et al. [73] 

1.3 

168 

17 

Simon et al. [72] 

1.43 

177 

18 

Simon et al. [72] 

1.53 

140 

14 

Simon et al. [72] 

1.75 

202 

40 

Simon et al. [72] 

2.34 

222 

7 

Delubac et al. [78], Ding et al. [79] 


1 


Table 2: Values of 

dA(z,) 

for different values of zbao- Here, Dy = 


1 

^ is the dilation 

Dv{zbao) 

ttAyZ) 

scale [86]. 



zbao 


dA[z*) 

Dv(zbAo) 


Reference 

0.106 


30.95 ± 1.46 

Beutler et al. [81] 

0.2 


17.55 ± 0.60 

Percival et al. [82] 

0.35 


10.11 ± 0.37 

Percival et al. [82] 

0.44 


8.44 ± 0.67 

Blake et al. [83] 

0.6 


6.69 ± 0.33 

Blake et al. [83] 

0.73 


5.45 ± 0.31 

Blake et al. [83] 
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Table 3: Best fit values of wq, oji and the minimum values of corresponding to the joint analysis 
of SN la + H(z) + BAO/CMB dataset with different choices of flmo- 


Name 

^mO 

Uq 

Ui 

Constraints on uq and ui 
(within Icr confidence level) 

Am^n 

GPL model 

0.26 

-1.07165 

0.424314 

-1.228< Wo < -0.9161, 
-0.3951< wi < 1.23 

579.536 


0.27 

-1.04934 

0.132091 

-1.22< Wo < -0.8766, 
-0.7968< wi < 1.062 

578.705 


0.28 

-1.02831 

0.103564 

-1.199< Wo < -0.8572, 
-0.8023< wi < 0.9991 

588.908 

JBP model 

0.26 

-1.07798 

0.6 

-1.301< Wo < -0.8525, 
-1.064< wi < 2.248 

579.667 


0.27 

-1.06072 

0.261937 

-1.293< Wo < -0.825, 
-1.518< wi < 2.003 

578.692 


0.28 

-1.05974 

0.024530 

-1.292< Wo < -0.8242, 
-1.789< wi < 1.737 

578.224 

BA model 

0.26 

-1.05344 

0.202546 

-1.182< Wo < -0.9232 
-0.2033< wi < 0.6094 

579.612 


0.27 

-1.04272 

0.059557 

-1.182< Wo < -0.9043, 
-0.4118< wi < 0.5235 

578.712 


0.28 

-1.06739 

0.018632 

-1.213< Wo < -0.9202, 

578.276 


-0.4808< wi < 0.515 


Table 4: Best fit values of Ag, flmo and the minimum values of corresponding to the joint 
analysis of SN la + H(z) + BAO/CMB dataset with different values of a. Note that u^{z = 0) = 
—As is the present day value of the EoS parameter for the GCG model. 


a 

A. 


Constraints on Ag and hlmo 
(within Icr conhdence level) 

\.min 

-0.5 

0.9 

0.240784 

0.8979< A, < 0.9462, 

0.20< fi^o <0.3396 

585.915 

-0.55 

0.9 

0.243458 

0.8976< A, < 0.9467, 

0.20< nmo <0.3432 

586.196 

-0.6 

0.9 

0.245692 

0.8997< A, < 0.9446, 

0.20< fl^o <0.3459 

586.344 










18 



(Oq 



a)Q 



COq 


Figure 7: Plot of la and 2a confidence contours on ujq — ui parameter space for the CPL 
parametrization (top panel), JBP parametrization (middle panel) and BA parametrization (bot¬ 
tom panel) respectively. In this graph, Xmin indicates the minimum value of corresponding to 
the best fit values of ojq and oji for the SN la + H(z) + BAO/CMB dataset, as indicated in the 
frames. The fixed value of fimo is also indicated in the frame. 

analysis is valid upto z > —1. The GCG behaves like dark energy = —0.9 > —1) at late time 
and its equation of state also settles to a value close to —1 in the far future. So, like the GGG 
model, the BA model will also avoid the finite-time future singularity [98, 99, 100, 101, 102, 103]. 
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A. 


Figure 8: Plot of la and 2a confidence contours on Ag — parameter space for the GCG 
parametrization. In this graph, Xmin indicates the minimum value of corresponding to the best 
fit values of Ag and for the combined dataset (SN la + H(z) + BAO/GMB), as indicated in 
the frame. This is for a = —0.5. 

It has also been noticed for the CPL and JBP models that the EoS parameter diverges in the 
hnite future {z ~ —1) as expected from theoretical expressions of u^p^z) for these models. 

In the limit, a(= oo, we have p</,(a) —)■ oo and p</,(a) —)■ oo for CPL as well as JBP 

models. This implies that the universe (according to the CPL and JBP models) will end up 
in Big Rip singularity [?], where the phantom energy density becomes very large in hnite time 
and overcomes the gravitational repulsion. Despite this future divergency problem (at z ~ —1), 
the CPL & JBP parametrizations have several advantages such as, they can probe the past 
evolutionary history of the universe, have a bounded and well behaved behavior for high redshifts, 
projects a simple two parameter (wq, a;i)-phase space with a good physical interpretation of ujq 
and Ui, etc. In addition to the simplicity, this kind of parametrizations (CPL & JBP) also have 
the advantage of giving hnite in the entire range —l<z< oo. Furthermore, we would like to 
mention here that the two-parameter models, adopted in this paper, provides a simple choice to 
obtain constraints from observational data. 


4 Conclusion 

It is well known that the late-time accelerating expansion of the universe can be described by a 
scalar held. For this reason, in this present work, we have discussed about four diherent types of 
non-canonical scalar held models with varying dark energy EoS for understanding the observed 
cosmic expansion. As a time-dependent EoS plays an important role for understanding the na¬ 
ture of DE, so we have considered four phenomenological parametrizations of dark energy EoS. 
The dynamical features of each models are analyzed, such as the evolutions of the deceleration 
parameter q{z) and the density parameters (D<^ and Dm)- The resulting cosmological behavior is 
found to be very interesting. 

For all the toy models, it has been found that the deceleration parameter q{z) indicates an 
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3 





Figure 9: The top and middle panel represents the plot of the EoS parameter vs. redshift using 
the best fit values of (uq, Ui) and Q^no = 0.27 (see Table 3). The bottom panel corresponds to the 
evolution of for the GCG model. This plot is for the best fit values of Ag and different values 
of a (see Table f); a = —0.5 (thick curve), a = —0.55 (dashed curve) and a = —0.6 (thin curve). 
The intersection point of the dotted line indicates uj^ = —1 at z = 0, the AGDM case. 


early deceleration followed by a late time acceleration of the universe (see figures la, 2a, 3a, 4a). 
We have also shown the evolution of density parameters and it is found that the results are in 
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good agreement with recent observations [64, 65]. 

We have also compared our theoretical models with the observational data coming out of 
the latest SN la, Hubble parameter, BAO and CMB measurements. For this purpose, we have 
written the Hubble parameter H{z) in terms of observable parameters {z, Hq & Hmo) and the 
corresponding model parameters for each DE parametrizations. We have obtained the best £t 
values of the parameters uq and ui by fixing the value of Hmo to 0.26, 0.27 and 0.28 (shown in 
Table 3). It may be important to mention here that the values of parameters of the model which 
were chosen for analytical results are well fitted in the Icr and 2cr confidence contours for each 
parametrizations. We have found that the oo^ = —1 crossing feature is also allowed by the SNIa 
+ H(z) + BAO/CMB dataset for the CPL, JBP and BA models with its present best-fit EoS 
parameter, wq < — 1 (as presented in Table 3). This is consistent with the results obtained by 
several authors [44, 93, 94, 95, 96, 97]. On the other hand, we have also found that the standard 
ACDM model is still compatible at the la confidence level for these toy models. It should be 
noted that the GCG model has been studied by many authors for the parameter range, 0 < a < 1 
[17, 18]. However, in this work, we have obtained the best fit values of the parameter Ag and 
Gmo by fixing the value of other parameter a within the range, —1 < a < 0 (shown in Table 4), 
as a < 0, is in good agreement with the work of Sen and Scherrer [104], Gong et ah [105] and 
Hazra et ah [41]. It has also been noticed from table 4 that the range of the allowed values of 
flmo match well with the previous results obtained by Riess et ah [106] and Sahni et al. [107]. 

However, as discussed in the previous section, the GPL and JBP models lose their prediction 
capability regarding the future evolution of the universe. We have also shown that it is very 
difficult to distinguish the GGG model from a AGDM in the near future and hence we need more 
investigations to constrain dark energy models more tightly. Obviously, we can not yet say which 
model is better as compared to other models by the analysis of the combined dataset (SN la -|- 
H(z) -|- BAO/GMB). We hope that the next generation observational data can provide more tight 
constraints on EoS parameter to enrich our understanding regarding the nature of dark energy. 
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